Efficient quantum memory and entanglement between light and an atòmic ensemble 

using magnètic fields 
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We present two protocols, one for the storage of light in an atòmic ensemble and the subsequent 
retrieval, and another one for the generation of entanglement between light and àtoms. They rely 
on two passes of a single pulse through the ensemble, Larmor precessing in an external field. Both 
protocols work deterministically and the relevant figures of mèrit - such as the fidelity or the EPR 
variance - scale exponentially in the coupling strength. We solve the corresponding Maxwell-Bloch 
equations describing the scattering process and determine the resulting input-output relations which 
only involve one relevant light mode that, in turn, can be easily accessed experimentally. 

PACS numbers: 03.67.Mn, 32.80.Qk 



INTRODUCTION 



Recent years have seen significant progress towards an 
efficient quantum interface between light pulscs carrying 
quantum information and atòmic ensembles suitable for 
storing and processing this information. Two approaches 
based on probabilistic photon detection or on de- 

terministic homodyne measurements 0, Q have been de- 
veloped. Of particular importance in the context of quan- 
tum information are means to swap the state of light and 
àtoms - enabling a quantum memory for light - and to 
crcatc Einstein-Podolsky-Rosen (EPR) type of entangle- 
ment of light and àtoms - the bàsic resource for quantum 
teleportation. 

Concerning the quest for a quantum memory, an im- 
portant experimental advancc was the recent demonstra- 
tion of the storage of weak coherent light pulses in àtoms 
P|, based on a Quantum Non Demolition (QND) interac- 
tion, measurement of light and feedback on àtoms. How- 
ever, reliable retrieval of the stored state by means of 
the same protocol would require the use of short pulses 
of squeezed light which are difficult to couple to atòmic 
ensembles in an efficient way. The design of less demand- 
ing protocols for storage and retrieval of states of light 
remained a challenge, also from a theoretical perspec- 
tive. Several protocols have been put forward, all relying 
on múltiple passes of light through the atòmic ensemble 
H BS ll E3, E3- Thc most efficient of these schemes, 
complying with the experimental requirement to use Lar- 
mor precessing atòmic spins, require eight passes of a sin- 
gle pulse 01 or two pulses each crossing twice an atòmic 
celi j 1 1 1 . In this paper we present a protocol, which con- 
sists of only two passes of a single pulse and achieves a 
state exchange of light and àtoms scaling exponentially 
in the coupling strength n, defined operationally as the 
signal to noise ratio of the undcrlying QND interaction. 
This scheme allows one to perform the complete transfer 
of a quantum state of light onto àtoms and back under 
modest experimental conditions, as we show for both, co- 
herent states as well as arbitrary superpositions of vac- 



uum and a single photon Fock state. 

Moreover, the same double pass setup serves with a 
slightly changed geometry as a deterministic source of 
EPR entanglement between light and àtoms. The en- 
tanglement scales thereby again exponentially in k. To- 
gether, these two protocols add to the growing toolbox for 
quantum information processing with room temperature 
atòmic vapors, which has already provided the possibility 
to entangle two atòmic ensembles via a Bell-measurement 
on two Larmor precessing spins [T3 | . In combination 
these tools undoubtedly pave the way towards numer- 
ous relevant applications, of which the demonstration of 
a complete quantum memory and quantum teleportation 
are just the most immediate. 

To be more specific, the setup of both protocols con- 
sists of an ensemble at room temperature in a cúbic glass 
celi. It is placed in an external magnètic field with large 
spin polarization along the axis of this field, such that 
the transverse spin components precess at frequency fi. 
A coherent pulse is directed through the atòmic samplc 
such, that it crosses it twice under an angle of 90 de- 
grees in the plane orthogonal to the axis of the magnètic 
field. The length d of the loop in the optical path is 
small, such that Larmor precession is frozen on a time 
scale d/c <§C f2 _1 , but the pulse length is large as com- 
pared to the Larmor period, T ^> Í2 . Under these 
conditions and the assumption that flT 3> k 2 , which is 
well fulfillcd in current experiments, we carefully solve 
the Maxwell-Bloch equations describing the dynamics of 
this scattering process. We identify the relevant light 
modes, which can be stored and retrieved or get entan- 
gled with àtoms and characterize their temporal profile. 
The central frequency of these modes lies at the upper 
or lower sideband of the carrier frequency, which is to be 
expected given the splitting of ground state levels of SI, 
and their slowly varying amplitude is exponential of thc 
form exp(±K 2 í/2T). The modes can thus be easily ac- 
cessed. Note that this setup is, apart from the magnètic 
field, similar to the one treated in 0] . It is precisely the 
presence of the magnètic field what enables us to achieve 
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FIG. 1: Relevant internal levels with quantization along x. 
Thick arrows represent the strong coherent field in x polariza- 
tion, thin arrows indicate the quantum field in y polarization. 



our results with the simple setup described above. 

The rest of the paper is organized as follows. In sec- 
tion [n] we introduce the bàsic idea of our protocol and 
summarize the central results. In section ITITl and IW1 we 
provide the detailed derivation for the quantum memory 
and the EPR source respectively. These sections are sup- 
plemented by two appendices. Finally, section [V] deals 
with sources of noise under realistic conditions. 



II. BASIC IDEA AND CENTRAL RESULTS 

In the following we consider a cúbic atòmic ensemblc 
at room temperature, which is placed in a magnètic field 
and interacts with a pulse of light propagating along z. 
The atòmic sample is assumed to be spin polarized along 
x, while the magnètic field is orientated along the op- 
posite direction. The pulse of light consists of a strong 
coherent i-polarized component of central frequency lüq, 
which is detuned by A from the atòmic transition, and 
a copropagating quantum field in y polarization. The 
beam's waist is assumed to cover most of the samples 
cross section. Àtoms have a relevant internal structure as 
shown in figure^ With x being the quantization axis, the 
classical light field drives the m — ±1/2 — * ml = ±1/2 
transitions, while the copropagating quantum field cou- 
ples to m = ±1/2 — > m' = ±1/2. In the case of a 
dominant ground state population of m — 1/2 levels, 
creation and annihilation operators of collective atòmic 
excitations can be defined by = £j| — 1/2)(1/2|/-\//Va 
and b, respectively, where Na is the total number of 
àtoms in the ensemble. Creation of an atòmic excita- 
tion will then be accompanied by the absorption (emis- 
sion) of a photon at frequency tt>o ± íí, (ujq — íl), that 
is, at the upper (lower) sideband, where Í7 is the Lar- 
mor frequency. Note that only the polarization, and not 
the energy of the sideband photons are relevant, so the 
notion of upper/lower sideband is rather arbitrary. Al- 




FIG. 2: Setups for having (a) a beam splitter or (b) a two 
mode squeezing like dynamics. (c) and (d) show the eflective 
transitions. 



though we will finally deal with light interacting with 
àtoms in free space, it is instructive to consider first 
the case, where àtoms are placed inside a cavity sup- 
porting both sideband modes. Related setups employ- 
ing cavities are considered in [ÏÜ. IÏ4| . We assume in 
the following that the cavity life time is much smallcr 
than the Larmor period SI -1 and let the creation oper- 
ators for the upper and lower sideband be given by 
and a, respectively. In the dispersive limit, the effec- 
tive Hamiltonian describing the interaction is given by 
H cx (bal is — b^aí ± h.c), where the signs follow from 
Clebsch-Gordan coefficients. Note that if ground state 
levels were degenerate, such that a us — a/ s = a, the 
Hamiltonian would be H cx (fe — b^)(a — a^), which is 
well known from the theory of quantum non-demolition 
(QND) measurements of atòmic spins. Including Zee- 
man splitting, the interaction consists of a passive and 
an active part, H = H pas — H act , where the passive 
part is a beam splitter Hamiltonian H pas cx b aj s + h.c. 
and acts only on the upper sideband, while the active 
part H act cx b^a, s + h.c. can be identified with a two- 
mode-squeezing interaction, which involves exclusively 
the lower sideband. Now, either of these two interactions 
can be selected in one of the setups shown in figure |2K 
or|2h- The interaction in every second pass will again be 
given by H but with phase changes a; s ( us ) — > ïO/ s ( us ) , due 
to the A/4 wave plate, and b — > ±ib, due to the change of 
the direction of light propagation, where the upper sign 
holds for setup in figure^ and the lower for[2). The re- 
sulting Hamiltonian is H' ± cx ±(60^ + £^a| s ± h.c). To- 
gether, we get for setup|21i an interaction H+H' + = H pas 
and for H + H'_ = —H act . In either setup one of the 
two A-type transitions in figure ^ is canceled by interfer- 
ence, and one is left with the transitions shown in figures 
|2t and 12b Note that these configurations remind of the 
Raman scattering processes put forward in |l5j| for the 
realization of a quantum repeater. Without a cavity, in 



3 



setups as shown in figuro [j3 the effects still persists, as 
we will show by solving the corresponding Maxwell-Bloch 
equations. In contrast to the dynamics inside a cavity, 
where Larmor precession of the atòmic spin is not crucial, 
it is well so for propagation in frec spacc. This can be 
understood by noting that both setups shown in figure[3 
possess a certain asymmetry in how the two transverse 
spin components in y and z direction are affected by light. 
This was calculated in great detail in 01 > where amongst 
others the setup of figure[21was examined without mag- 
nètic field. We emphasize that Larmor precession helps 
to remove this asymmetry. 

In the rest of this section we collect the results for 
both, the quantum memory and the two mode squeezing 
protocol. This will be done in the language of canonical 
operators xa = (b+ w)/y2 and pa — — tf)/y/2 and 
likewise for light, since solutions to Maxwell Bloch equa- 
tions are more conveniently derived in this formalism. 

Quantum memory Within the memory scheme, figure 
the transfer of a quantum state of light onto àtoms 
or vice versa approaches perfect mapping exponentially 
in the coupling strength. We have 



(a) 



X A 
out 

Pa 



= P. 2 



Pa 



in 

Pl+ 



r . n . . * i in/ out n in/ out 

tor tnc write-m procedure, wnere x A and p A are 
the atòmic input / output quadratures of the scheme and 

in/out i in/out r , n . n . n . n T , 

x L + and p L ^_ refer to the write-m light mode. It 
lies at the upper sideband (according to the configuration 
considered above) and is modulated by a slowly varying 
envelope with an exponential profile, which is a propa- 
gation effect. For the retrieval the inverse accented light 

n 'in/out n An/out . i i i 

mode x L _ and p L _ is used and we have 



x 



in 

Pa 



'in 

Pl+ 



where accents indicate quadratures referring to the read- 
out pulse. Note that for large n the state exchange is 
perfect. It is remar kable that both pairs of input-output 
relations have a form which reminds of a decoherence pro- 
cess, with the important difference that we have modes 
in place of Langevin noise operators, which can be con- 
trolled at will. The fidelity for the complete state transfer 
- write in and subsequent retrieval of a state of light- is 
given in figure^Ja) and (b) for coherent input states and 
light qubits respectively. 

EPR source The active version of the protocol, fig- 
uro [3Jd, generates correlations between àtoms and light, 
which grow exponentially in the coupling. One can define 
interspecies EPR modes 

1 / - N 1 / - N 

Xl = —^(XA~ PL+j, pi = -j=(jpA + XL + ), 

1 / - N 1 / - X 

X2 = —^(XA+PL+j, p 2 = —^(PA- X L + ). 

where ïj j+ , p^ + refer to a light mode, which resembles 
the mode xl+, Pl+ introduced above apart from the fact 
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FIG. 3: Schemes for realization of a quantum memory (a) 
and a source of EPR entanglement (b) . (a) In the first pass 
a pLPA-interaction oceurs. Subsequently the pulse is sent 



through a 



plate, which interchanges xl and pl- The 



pulse is reflected back onto the sample. This happens at a 
timescale much shorter than the Larmor precession of the 
àtoms. Therefore the transverse components of the collec- 
tive spin can be assumed to remain in their place to a very 
good approximation. Finally the pulse passes the àtoms along 
y. Due to the changed geometry the atòmic quadratures are 
also interchanged. pa — * xa and xa — * —Pa, which means, 
that the light field couples to xa in it's second passage, hence 
leading to a xlxa -interaction. In (b) the changed geometry 
introduces a different sign in the exchange of atòmic quadra- 
tures, which leads to a —xlxa interaction in the second pass. 



that the lower sideband is involved instead of the upper 
one. x\ and pi are squeezed, while X2 and p\ are anti- 
squeezed, 



(Axi) 
(Api) 



(A P2 ) 2 
(Ax 2 ) 2 



—i k 

where z = cosh (e~). The EPR variance of the gener- 
ated state is depicted in figure 

The results presented above will be derived in the fol- 
lowing sections. We remark that each protocol can be 
realized involving either the upper or the lower sideband. 
The sideband mode involved can be changed by either in- 
verting the ground state polarization or changing the ori- 
entation of the magnètic field. Losses will be considered 
in section Ivl and it will be shown that the proposed pro- 
tocols are robust against the dominant sources of noise. 
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III. QUANTUM MEMORY 

The following calculation will be done by means of 
canonical operators. For àtoms canonical variables are 
defined by means of the Holstein-Primakoff transforma- 
tion and approximation |l6j | . Via the Holstein-Primakoff 
transformation spin-eigenstates are mapped onto har- 
mònic oscillator-eigenstates. The initial coherent atòmic 
spin state is treated as harmònic oscillator ground state. 
The Holstein-Primakoff approximation allows one to de- 
fine the canonical atòmic variables xa and pa corre- 
sponding to the y and z component of the collective angu- 
lar momentum J, xa = J y / \/ {J x ) and pa — J z / y/ (Jx)- 
The light field in y-polarization is described by spatially 
localized modes 



X L 



(r) = -J= / M^e-^°-^ r/c + h.c), 

V47T J b 



du^e-^ -^ - h.c), (1) 



where the range of integration b is a small bandwidth 
around the carrier frequency ujq containing Í7. The spa- 
tial argument r refers to the distance along the optical 
path shown in figure |3| and we have [xl{t),pl{t')\ — 
icS(r — r'), where c is the speed of light and the width 
of the delta function is on the order of c/b Within this 
description the Hamiltonian for the off-resonant scatter- 
ing interaction takes the form H oc p wa |l7j |. Detailed 
descriptions can be found in and . 



A. Write-in 

The double-pass interaction in setup can be de- 
scribed by 



H = H a 



Hi 



light 



H àtoms = —{x 2 a + Pa) refers to Zeeman-splitting of 
the atòmic ground state causing Larmor precession of 
the transverse spin components represented by xa and 
PA- The interaction terms V\ and V% account for the 
off-resonant scattering interaction in the first and second 
passage of the pulse respectively. They are given by 



V, 



Hk 



PaPl(0) and V 2 



Kk 



x A x L (d) 



where T is the duration of the pulse. V\ was already 
introduced. V2 basically describes the same kind of in- 
teraction, but due to the changed geometry in the sec- 
ond pass atòmic quadratures are interchanged pa — > xa- 
Since the beam is sent through a quarter wave plate 
between it's passes through the atòmic sample, light 
quadratures are interchanged as well pl — > xl- The 
arguments of the light-operators in V± and Vi indicate 
that the first scattering interaction occurs at r = 0, 
while the second interaction happens after the light has 



travelled some distance d in the small loop between the 
mirrors. The length of the làser pulse is hereby sup- 
posed to be large compared with the distance within 
the loop. In typical experiments pulses of a length of 
several hundred km are used, therefore the pulse en- 
counters itself in the sample 27]. Hu g ht represents free 
propagation of light. It acts on light quadratures like 
d t XL(r) = j-[H lMjht ,x L (r)] = ~cd r x L (r), which is a suit- 
able approximation for the light modes defined in ft}. 
Evaluating the Heisenberg equations gives 



d t x A (t) 
dtPA(t) 
(d t + cd r )x L (r, t) 
(d t + cd r )p L (r, t) 



PL(0,t) , 



Íïpa(í) 



-üx A (t) =x L (d,í) 

vi 



T PA(t)S(r) , 
~j=XA(t)8{r - d) 



By performing the variable transformation £ = ct — r wc 
obtain the Maxwell-Bloch equations 



d t x A (t) 
d t PA{t) 

dtPLÍU) 



-Qxa(í) -j=x L (ct - d, í), 



PA (t)S(ct - o, 



--j=x A {t)5{ct-Í-d). 



(2) 
(3) 
(4) 
(5) 



Light modes in new variables are denoted by a bar 
^i(^jí) — xl{cí — £, í). The light variable argument Ç 
refers to a coordinate system which is fixed on the light 
pulse. It allows one to denote easily particular pieces on 
the pulse. At a certain instant of time £ labels the pieces 
according to their position starting with the piece, which 
enters the atòmic sample first. 

This set of coupled diffcrcntial equations has now to be 
solved. As a first step we treat the equations for light. 
In the first pass a PlPa -interaction occurs and xl picks 
up some pa contribution. The delta function in (0} re- 
flects the fact that a certain piece £ of the pulse gets a 
contribution from the atòmic state at t = £/c (which is 
the instant of time the piece in consideration passes by). 
In the second pass a xlxa -interaction occurs, and the 
atòmic x -quadrature is written onto pl- A piece £ of 
the pulse, which interacted with p^ at time £/c gets a 
contribution from x l after it has traveled a distance d in 
the loop. Therefore the atòmic x quadrature is picked up 
at t = £/c + d/c which is indicated by the delta- function 
in equation (JSJ) ■ By integrating equations (@J and JHJ| for- 
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mally these delta functions turn into Heaviside functions, 

x L {t t) = x L {Ü, 0) / drp A (T)S(cr - Ç) 

Vi Jo 

= afc(£,0)+-!LpA(É/c)e(t-e/c), 

V -í 



Pi (6 *) = 0) - ^= / ^ a (t)<5(ct - £ - d) 
vi io 



xa{Ç/ c + d/c)0(í - - d/c). 



Now p L (ct, t) and x_l (cí — d, í) are calculated, since these 
expressions have to be substituted into the atòmic dif- 
ferential equations @ and (J3J. The fact that the argu- 
ments are difïerent for Xl and pi can be understood by 
considering the processes going on in the course of the 
double pass scheme. During the pip^-interaction in the 
flrst passage xa picks up some pl contribution. If wc 
consider this process at a certain instant of time í, the 
relevant piece of the pulse is the one passing r = 0. It is 
denoted by £ = cí — r = ct. pa is acted upon by xl in the 
second pass by the piece of the pulse which passes r = d 
at time í. So it gets a contribution from p^{ct — d, i). 
One finds 



Xl(cí — d, í) 



p L (ct,t) 



x L (ct-d,0) + —i=Pa (t - d/c) 9 (d/c) 



x L (ct - d, 0) + — p A (t - d/c), 
vT 



p L {ct, 0) - -=x A (t + d/c)Q(-d/c) 
VT 

p L (ct,0). 



Note that p L (ct,t) is conserved. This feature is due to 
the time-delay in the loop and will turn out to be crucial 
for the characteristic exponential behavior of the wholc 
scheme. After inserting these results into © and © the 
atòmic differential equations read 



a t x A (t) = sip A {t)+-=js L (ct ) o), 

vi 



d t p A {t) = -Qx A (t) ^=x L (ct - d,0)-—p A {t - d/c). 

Now we assume d/c <C Í2 _1 , such that the elapsed time 
during the run in the loop is definitely much shorter than 
any other relevant process. d/c can be assumed to be of 
the order of ns while àtoms rotate slowly with a Larmor 
period of the order of fis. With this approximation 



dt [pA(t) 



= < n 



o i 

-i o 



K 



o o 

1 



( x A {t) 

\PA(t) 



Pi(cí,0) 

y/T \—XL(ct,0) 



(6) 



This differential equation consists of a homogeneous part 
and a driving term. The first term of the homogeneous 



part - being proportional to the Larmor frequency - re- 
flects the fact that àtoms turn with Í7 in the external 
magnet field. The second term in the homogeneous part 
represents damping of pa ■ Although only one quadrature 
is damped, the effect is distributed among both quadra- 
tures by Larmor precession. This leads to a symmetry 
between x and p, which is a characteristic feature of our 
proposal. The solution to the differential equation is 



x A (t) 

PA(t) 



A(t) 



+A(t) 



x A {0) 
Pa(0) 



where A(t) = e 



Gi 



G = n 



dT A~ 1 {t) 



1 

-1 



Pl(ct,0) 
-x l {ct,0) 



00 
r i 1 



is the homogeneous solution. We supposc flT ^> k 2 , 
which matches experimental conditions, since typically 
Í7T rj 300 while n 2 is of order unity. With this 
assumption 

A(t) = e-&R-\t) , 

where iï _1 (í) is an orthogonal matrix, 

fí -ln\ _ ( cos(fií) sin(Oí) \ 
( ) ~ \ - sin(Oí) cos(Oí) ) ■ 

The inverse is taken for later convenience. Therefore the 
atòmic time evolution is given by 



x A (t) 
PA(t) 



in 
X A 

in 

Pa 



(7) 



+e~^rR~ 1 (t)— = / dre^R(T 



Pl(ct,0) 
-xl(ct, 0) 



Now the atòmic output quadratures x° A ut — x A (T) and 
pout _ p A (T) can be directly written down. With the 
assumption QT = 2itn for some natural number n, 



x A 

out 

Pa 



— BZ. / X 



A 

in 

Pa 



+er^—= / dt e^R{t) 



Pi(cí,0) 
-x L (ct, 0) 



The atòmic output-quadratures consist of some atòmic 
input contribution which is damped exponentially with 
k 2 and an additional light contribution which they pick 
up during the scattering interaction. The definition of 
the appropriate light-mode can be taken from this result 
right away, 



■ L· l+ 

in 

Pl+ 



Tvc 



dt e%rR{t) 



Pl(cí,0) 
-x L (ct,0) 



(8) 
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where the prefactor assures normalization such that 
[x™ + ,p™ + ] = i. This new defined light mode is essen- 
tially the upper sideband mode x us , p us , which is given 
by 



quadratures after the whole pulse run through the atòmic 
sample 



in 
fus 



1 



—= dt R(t) 



-x L (ct,0) 



(9) 



The only difference is given by the fact that the stored 
mode x™ + , p™ + is defined with a slowly varying envelope 
of the form exp(+K 2 í/2T). The index " + " refers to the 
sign of the argument in this exponential function (later on 
we will also have to deal with corresponding " — " modes) . 
With use of |[SJ the atòmic input-output relations can be 
written in a compact form 



~,out 

X A 

out 

Pa 



e 2 



in 

Pa 



VT 



■ L· L+ 
in 

Pl+ 



(10) 



These equations describe the write-in process for a sig- 
nal, which is encoded at the mode described above. Re- 
markably, mapping of such a quantum state of light onto 
àtoms approaches perfect read-in exponentially in the 
coupling strength. This arises from the fact, that in the 
coursc of the double pass scattering interaction xl picks 
up some contribution from the atòmic p -quadrature, 
while pl in contrast is conserved. Therefore we do not 
get a rotating term in the bàsic differential equation (jüj , 
which would lead to sines and cosines in the solution, 
as we would expect for a beam splitter like interaction, 
but an exponential effect, which is characteristic for the 
setup. 



To evaluate this expression in terms of input-operators 
the integrated versions of equations (@J and © 

MiJ) = é L (£,o) + -^4(£/ c )e(í-£/ c ), 



T 

Vdtt) = J> L (Ç,0)--±=x A (Z/c)Q(t-t/c). 
are used. Therefore 



'out 

Pl- 



VTVl - e-* 2 
%(cí,0) 
-xl (ct, 0) 



dt e~^rR{t) 



Now the atòmic time evolution iQ has to be inserted. 
The resulting expression can be simplified by inter- 
changing the order of the double-integral J Q T dt J * dr — > 

lo d T f T dt. With help of equation © the read-out out- 
put can then be written as a sum of an atòmic contribu- 
tion and some contribution from the plus-mode. 



'out 
X L- 
'out 

Pl- 



Pa 



PÏ+ 



(11) 



B. Read-out 

In order to perform the read-out, a pulse of light has 
to be sent through the double-pass setup, just like for 
the write-in procedure, but since we are now looking at 
the reverse process, the appropriate light mode for this 
task has to be accented in an inverse fashion. While 
in the write-in process the rear part of the pulse was 
emphasized, now the front part of the pulse has to be 
weighted in order to pick up atòmic information best. 
As the exponent in the mode definition is negative this 
read-out mode will be denoted by a minus sign. Since we 
now deal with a new beam of light which is independent 
from the write-in pulse, read-out beam variables carry an 
accent, 



Pl- 



VtVT 



dt e 



l R(t) 



%(cí,0) 
-x L (ct,0) 



with a new normalization constant k/(VtVÏ — e~ K - À ). 
The input-output relations for this mode can be derived 
by changing the time argument of light operators from 
to T, reflecting the fact that we now look at the light 



Note that this expression resembles the formula for the 
write-in procedure with the roles of light- and atòmic 
modes interchanged. 



C. Fidelity for the complete state transfer 

The fidelity for the complete state transfer is given by 
the overlap of the initial input state and the final output 
state after storage and subsequent retrieval. By inserting 
the output of the write-in procedure l|10(l into the read- 
out equation Ullfl one obtains 



fin 
X L _ = 

fin 

Pl- = 



-(1 
(1 



L + 

Pl+ 



■VT 



e-^x'X 



■VT 



e- 2 p% 



- e 
e~ 



x L-i 

pT- 

(12) 



For infinite coupling k 2 the original input-signal is re- 
trieved within the final quadratures of the read-out pulse 



L· L + 



and p L ™i — 



while the noise terms 
(atòmic and read-out beam input contributions) vanish. 
The quantum state to be stored is supposed to be un- 
known. It is assumed to be taken from a certain set of 
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possible input-states. In the following two subsections 
we will consider coherent input states and light qubits 
respectively. We will first calculate the fidelity for a sin- 
gle state transfer and take the average over the complete 
set of possible input states in the next step in each case. 
The results will be compared to the corresponding clas- 
sical limits, i.e. the maximum average fidelity, that can 
be achieved by classical means [2(J, W& Hj, |23|, UÈ liÉ, ■ 



1. Fidelity for coherent input states 

We first consider storage of a coherent state of light. 
The overlap between an initial state with quadratures 
x l+ > P % l+ an d the final state with é , p^"™ is given by 



F, 



coh 



[l + 2(Afi m )2][l + 2(A^ 

(<4+>+<4'"» 2 «pí+>+^{'">) 2 



l + 2(Ax{ i _™)2 l + 2(Ap{^™)2 



(13) 



The expectation vàlues and variances of the final light 
state follow directly from (|12[) . Since the àtoms and the 
read-out plus mode are initially in a vacuum state we 
have (x{ m ) = -(1 - e- K2 ){xf + ) and (p£?) = -(1 - 
e^ K2 ){p l L n + ) , while the variances are given by (Ax^) 2 — 

(Ap£ m ) 2 = i, as one expects for a passive transforma- 
tion. Therefore 



where = -^{ x T+~ vT+) 1S t üe creation operator for 
a photon in the write-in mode. The write-in and read-out 
procedure is given by a passive transformation U 



(a ■ 
ia- 
ia- 



(3U a L ac) 
f3Ua[™tf)\v, 
f3a L l_ in )\vac) 



(14) 



where U\ vac) 



\vac) was uscd. Hcre a 



t fin 



fin 



W'l-) is the creation operator after mapping and subse- 
quent retrieval. It can be directly calculated, since the 
completo input-output relations for the light-quadratures 
are known. With use of equations (|12f) one finds 



a 



t fin_ 



where a^"\ 



f in 



(15) 

and à L ™ refer to the light state to be 
stored, the àtoms and the read-out mode respectively. 
The fidelity is given by the state overlap between \^f m ) 
and the optimal final state \^{,pt) = (a — (3a L ™)\vac). 
By inserting (|15fl into expression l|14|) Fq^n can easily 
be determined. One obtains 



fin\\2 _ 
opt / I — 



,| 2 + {l-e-« }\(3\ 2 )f 



The average fidelity is calculated by setting a = cos(|) 
and (3 — sin(|)e^ and integrating over the whole Bloch- 
sphere, 

F qubl tÍK) = -!- / <!(-) 



dej) sm(9)F qubit (6», 0) 



Now the average fidelity is computed by averaging over 
the complete set of all possible coherent input states. For 
this purpose the amplitudes x l £ + and p l £ + are assumed to 
be taken according to a Gaussian distribution centered 
at zero with a certain width n. 



1 



Figure 2t> shows this result. The maximal average fi- 
delity that can be achieved for qubit states by a classical 

J| is indicated by a cross. 



strategy Ff uUt 



F C ohin, k) 



F coh ((xt + ),(pf + ), K ) , 
1 

1 + e~ 2Kj2 n 



Figure 0Ji shows the average fidelity for different widths 
corresponding to mean photon numbers of the distri- 
bution. The corresponding classical limit Ff oh — 



2fj,|21j is marked by a cross on each curve. 



2. Fidelity for light qubits 



IV. EXPONENTIAL TWO MODE SQUEEZING 

The interaction which governs the squeczing scheme 
pictured in figure EHb), is given by 



H = H„ 



Hi 



light 



v 1 -v 2 



This Hamiltonian diffcrs from the one used in the mem- 
ory section just by a sign in the interaction term referring 
to the second passage. The pulse runs along —y in the 
second pass of the squeezing scheme (instead of y in the 
previous case) and sees therefore —xa- Hence we have 
the minus sign in front of Va for the new setup. 



Now the fidelity for light-qubits is calculated. The 
light-qubit input state is represented by 

l* in > = {a + f}a\£)\vac), 



A. Input-output relations 

The atòmic input-output relations can now be derived 
in complete analogy to section IÏII Al By evaluating the 
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FIG. 4: Average fidelity for write-in and subsequent read- 
out of a light state versus coupling k 2 . Crosses indicate the 
classical limit in each case. (a) Average fidelity for coher- 
ent light states according to distributions with different mean 
photon numbers (sòlid line: n—4, dashed line: n=8, dotted 
line: n=20) (b) Average fidelity for light qubits. 



Heisenberg equations as above we get 
x A {t)^ ( < " " -'- 



d, 



p A (t) 



= m 



o i 

-i o 



KT_ I 

T VO 1 



x A {t) 

PA{t) 



^/f \x L (ct,0) 



With the usual approximation k 2 <C 2ÍIT we obtain 



x A (t) 

PA{t) 



e^ r R~ í {t) 



'< A 

p'a 



-e^rü _1 (t) 



t 

dre~ 



X L (cT,0) 



These equations are in a significant way different from 
the atòmic time evolution Q in the memory scheme. 
Note first the signs in the arguments of the exponential 
functions. We now have exponential enhancement of the 
atòmic input instead of exponential damping. Further- 
more light is involved in form of a minus mode in the 
atòmic input-output relations because of the minus sign 
in the exponent within the integral. Note second, that 
the minus sign, which was present in front of Sl(cí, 0) in 
the memory scheme, does not appear in this case. There- 
fore the lower sideband 



in 
Pis 

in 
X ls 



1 



T , 



(16) 



is involved instead of the upper one (JSJ . Hence the minus 
mode showing up in the atòmic time evolution is defined 



slightly differently from the memory section 



Pl- 



vWÏ - er 



dte~^rR{t) 



P_ L ipt,Q) 
x L (ct, 0) 



With use of this definition and the assumption Í7T = 
2irn for some natural number n, the atòmic input-output 
relations read 



ryOUt 

x A 

OUt 

Pa 



Pa 



(17) 



Light input-output relations for this process can be de- 
rived in analogy to the procedure in section Ull Bl The 
inverse accented counter-part of the light mode used in 
the atòmic evolution is given by 



VTVe^ - 1 Jo dte2TR{t) \x L {cm 



We refer now to the same pulse as in JT7J, while in the 
derivation of the input-output relations for àtoms and 
light in section lïïïl two independent beams were consid- 
ered. We obtain 



~out I — V < 



1 



Pa 



4 pï- 



(18) 



Please note that the input-output relations l|17|) and i|18|) 
are active versions of (110(1 and pi(l respectively. 



B. Creation of entanglement 

As can be seen from the input-output relations for 
àtoms and light given in equations (|17fl and l|18|l respec- 
tively, correlations between àtoms and light are created 
which grow exponentially in the coupling strength. We 
define new modes appropriate to the type of correlations 
produced in the system by setting 

1 / - ^ 1 / - x 

Xl = -j={X A - PL+), Pl = -y={p A +X L+ ), 

1 / - N 1 / - > 

X2 = — (X A +p L+ ) : p 2 = —7=(PA ~ X L+ ). 



V2' 



V2' 



The corresponding variances can be calculated easily 
from JEJ and ÏÏB- We S et 



(A Xl ) 2 = (Ap 2 ) 2 = (Ve- 2 - 1 - 



.2 N 2 



2 = 



(A Pl ) 2 = (Ax 2 ) 2 = (Ve K2 + = 



2 s 2 



— 1 K 

with z — cosh (e~). We get a two mode squeezed 
state where x\ and p 2 are squeezed, while p\ and x 2 are 
antisqucezcd. In the limit of infinite coupling the state 
becomes an EPR state in which x A , p^ + and p A , Xl+ are 
perfectly correlated. For the state under consideration, 
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FIG. 5: Entanglement produced in the two mode squeezing 
scheme versus coupling k 2 . The entanglement is hereby mea- 
sured by the EPR variance Aepr- 




FIG. 6: Atòmic squeezing (Ap a ) 2 in db versus coupling re 2 
for optimal gainfactor g op t- The inset shows how the optimal 
gainfactor depends on the coupling. 



holds in the ensemble average, as is shown in 0, IÏ8j . 
With help of the atòmic input-output relations fïTf and 
the expression for the measurement outcome above one 
finds 



fb 

Pa : 



/e^ - 5 Ve K2 -ljp%+ (Ve K * -l-ge 



Thus the variance of this quadrature is given by 

2-, / „\ 2-, 



(Ap{ h f= [e^-g 



lV e» 2 - 1^ ~ + (V e« 2 - 1 - ge£\ ~. 



(Ap A b ) 2 is now optimized with respect to the gainfactor 
g. We obtain 

gopt 



e 2 



Ve* 2 - 1 + e K VI - e~ K2 



2e K " - 1 



(A^ b opt ) 2 = 



1 1 



2 2e« 2 - 1 



Note that the àtoms are left in a minimum uncertainty 
state, since 



1 



(Ax A y = ~(2e" - 1) = - 



1 1 



4 (Ar, f . b )2 



The amount of squeezing depending on the coupling n 2 
is shown in figureEJ 



the EPR- variance Arpb = \{Ax\ + Ap 2 ) — e~ 2z is an 
entanglement measure [26| . For separable states A epr — 
1. For inseparable states Aepr decreases with increasing 
entanglement. The amount of entanglement created in 
the scheme is shown in figureEl 



C. Spin squeezing 

The correlations created in the proposed scheme can be 
used to produce atòmic squeezing. This can be achieved 
by performing a measurement on the plus light mode 
and subsequent feedback onto the atòmic spin based on 
the measurement outcome. The squeezing protocol is 
symmetric with respect to the interchange of {xa, Pl} 
and {pa, Here squeezing of (Apa) 2 is illustrated. 

In order to aequire information about pa, xl+ has to be 
measured. The outcome of this measurement is governed 
by the operator equation 

x°£l = Ve* 2 - 1 p A l + e^xT-- 

If the measurement outcome is obtained p° A ut is dis- 
placed by an amount g where g € M is some gain 
factor. For this feedback procedure the operator identity 

r/ b — n out — n T- out 

Va — Va 9 x l+ 



V. CONSIDERATION OF NOISE 

We consider losses for both components of the protocol 

- àtoms and light - and treat them perturbativly within 
the Gaussian formalism. Concerning the atòmic sample 
we take transverse decoherence of the atòmic spin state 
at a rate of ^ into account. As in experiments atòmic 
vapor is usually contained within a glass celi, the domi- 
nant source of noise concerning light are reflection losses. 
These affect both, quantum variables and classical field 
and will be characterized by the reflection coefRcient r. 

A. Quantum memory with noise 

In this section we sum up results for write-in and read- 
out in the presence of losses. A detailed description of 
the generalized quantum memory scheme including noise 
is given in appendix El Consideration of losses leads to 
a modification of the original write-in mode. The gener- 
alized write-in quadratures preferred by the system are 
given by 

where w = rj/T + k 2 (1 — 2r)/T. Both sources of noise 

- reflection losses and spontaneous decay as well - give 
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« 2 K 2 



FIG. 7: Average fidelity with losses versus coupling. The 
atòmic decay rate rj and the reflection coefficient r both have 
a value of 7.5%. Crosses mark the corresponding classical 
limits. (a) Fidelity for coherent input states according to 
distributions with different mean photon numbers. (sòlid line: 
n = 4, dashed line: n = 8, dotted line: n = 20) (b) Fidelity 
for light qubits. 
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r r 

FIG. 8: Maximal attainable average fidelity for coherent input 
states according to a distribution with mean photon number 
n = 8 (a) and light qubits (b) versus reflection coefficient r 
for different atòmic decay parameters rj. (sòlid lines: rj = 5%, 
dashed lines: rj = 10%, dotted lines: r\ — 25%) The dash- 
dotted line and the cross indicate the classical limits. 



rise to a generalized exponent in the exponential modu- 
lation function. In addition to the changed envelope, the 
light mode appearing in the atòmic input-output relation 
is further disturbed: it lies no longer exactly at the up- 
per sideband, but contains a small contribution from the 
lower one, as can be seen by comparing the expression 
above to JÏÏJ and l|l(j|) . Since it is experimentally advan- 
tageous to encode the input signal at sideband modes, 
we define a generalized write-in mode (denoted by capi- 
tal letters) 

which takes full account of noise concerning the expo- 
nential modulation, but lies precisely at the upper side- 
band. (i.e. the small orthogonal contribution from the 
lower one is treated as noise.) To perform the read-out, 
the inverse accented counter part of this mode is mea- 
sured. The calculation of the fidelity is given in appendix 
IBI Figure [7| shows the average fidelity for write-in and 
subsequent retrieval versus coupling for r = rj = 7.5%. 
Plots (a) and (b) refer to coherent input states and light 
qubits respectively. The corresponding classical limits 
are marked by crosses. As illustrated by these graphs 
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r 

FIG. 9: (a) EPR variance Aepr versus coupling k 2 in the 
presence of losses. The reflection coefficient r and the atòmic 
decay rate are both chosen to have a value of 10%. 
(b) Optimized EPR variance versus coupling for different 
atòmic decay parameters. (sòlid line: rj — 5%, dashed line: 
rj — 10%, dotted line: rj = 25%) The inset shows how the 
optimal coupling k opt varies with r. 



losses decrease not only the quality of the state transfer 
for a given coupling strength, but limit also the attain- 
able fidelity. The crucial limiting factor in this scheme are 
reflection losses. Figure [S] shows the maximum average 
fidelity versus r for different vàlues of the atòmic decay 
parameter rj. Plot (a) shows results for coherent inputs, 
while plot (b) depicts the maximal attainable fidelity for 
qubits. The dash-dotted line and the cross indicate the 
classical limits in each case. Within moderate couplings 
fidelities well above the classical limit can be achieved, 
showing that the protocol is robust against the dominant 
sources of noise. 



B. Two mode squeezing with noise 

Consideration of noise within the two mode squeez- 
ing protocol is done along the same lines outlined in the 
section above. The entanglement created by the scheme 
in the presence of losses is depicted in figure HJa) for 
r = r] = 0.1. The EPR variance increases for higher 
vàlues of k 2 . An optimal value K 2 tpt exists for which 
the proposed protocol works best and a maximal amount 
of entanglement is generated. Figure EJb) shows the k- 
optimized EPR variance versus r, while the dependcncc 
of K op t on the reflection coefficient is given within the in- 
set. As can be seen from these plots atòmic decay plays 
a minor role. 

Spin squeezing can be performed with a lower and lim- 
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FIG. 10: (a) Spin squeezing in db versus coupling k 2 in the 
presence of losses. The reflection coefficient r and the atòmic 
decay rate both have a value of 10%. The inset shows how 
the optimal gainfactor g op t depends on the coupling. 
(b) Maximal spin squeezing versus reflection coefficient r for 
different atòmic decay parameters. (sòlid line: r\ — 5%, 
dashed line: rj = 10%, dotted line: r\ = 25%). 



ited quality in the presence of losses. In contrast to the 
ideal case the optimal gainfactor does not approach unity 
with increasing coupling but converges towards a higher 
value which depends on the amount of losses impairing 
the system. Figure lTUT a) shows the squeezed atòmic vari- 
ance in db and the dependence of g op t on the coupling for 
r = T) = 0.1. The maximal attainable squeezing ver- 
sus r is given in figure HUf b) for different atòmic decay 
parameters. 



VI. CONCLUSIONS 



In conclusion we propose two protocols based on a 
double-pass scheme for a single atòmic ensemble in a 
magnètic field. The first protocol provides an exponential 
scaling interspecies beam-splitter interaction. Therefore 
it is suitable for high hdclity storage and retrieval of an 
unknown quantum state under modest experimental con- 
ditions, as was shown for coherent input states and light 
qubits as well. The second protocol generates determinis- 
tically EPR entanglement between àtoms and light. The 
proposed protocols provide therefore the ingredients to 
realize a variety of interesting quantum communication 
protocols. They are also shown to remain experimentally 
feasible under realistic conditions. 
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APPENDIX A: CONSIDERATION OF NOISE IN 
THE MEMORY PROTOCOL 

Write-in 

Atòmic noise can be incorporated into the framework 
of section lïïïl bv including decay terms in Bloch equations 
© and © 



d t x A (t) = Qp A (t) + -j=p L (ct,t) - 7^,X A (t) 
d t p A (t) = -Qx A (t) ^=x L (ct- d,t) 



2T 



PA{t) 



fpA(t), 



rj 

jífx A {t), 



(Al) 



where f xA and f pA are Langevin noise operators with 
zero mean and (f{t)f{t')) = S(t -t')\. 

Each time light crosses one of the celi walls, reflec- 
tion losses occur. This happens four times. In the 
following we will consider coherent input states. In this 
case losses due to the very first crossing can be neglected, 
since these can be compensated by using a more intense 
input signal. In case of light qubit input states losses 
due to the first reflection have to be considered, which 
makes the resulting equations slightly more complicated. 
Since the derivation is analogous except for this point 
it won't be sown explicitly. Losses due to the second 
and third trànsit of a celi wall affect only the second 
scattering interaction. We take this into account by 
modifying the undisturbed equations for the light field 
quadratures to be inserted into IjAlll 



XL{ct—d,t) = xl{cí— d, 0) H — j=p A (t 
p L (ct,t) = p L (ct,0), 



by introducing light quadrature damping with a factor 2r 
(the factor 2 reflects the fact that crossing of a celi wall 
happens twice) and corresponding light-Langevin opera- 
tors f X L and f P L, and obtain 

x L {ct-d,t) = VI - 2r(x L (ct - d, 0) + A=p A (t - -j) 



p L (ct,t) = p L (ct,0). 



p L (ct, t) remains unchanged, since this quadrature affects 
the àtoms only in the first passage (during the plPa - 
interaction), which means that each part of the pulse con- 
tributes before it is subjected to reflection losses. There- 
fore pl is conserved as in the undisturbed case. The 
classical light field is impaired by reflection losses as well. 
Since the coupling strength of the scattering interaction 
is proportional to the amplitude of the classical field we 
have a reduced coupling k = \/í — 2r k for the second 
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(xlxa -) interaction due to the light crossing two celi 
walls before it's second passage. By considering this and 
inserting the expressions above into equations IjAljl we 
obtain 



d t XA (t) = n PA (t) + -^=p L (ct,0)-^X A (i) - 



fxA(t), 



dtPA(t)=-Qx A (t)—= y/ï^2Ï(x L {Ç,0) + ^ PA (t) 

Vi L V Vi 



2rf xL (t) 



2T 



PA{t)- 



f P A(t). 



We can ignore reflection losses arising in the very last 
trànsit through a celi wall, since the light field of the 
write-in beam is of no relevance after the second scatter- 
ing interaction. By neglecting the time delay d/c as in 
section lïïïl the atòmic differential equations generalize to 




1 
1 



k 2 {ï 



T 



p L (ct,0) 
-(1 - 2r)x L (ct,0) 




We introduce the abbreviation w = rj/T + k 2 (ï — 2r) /T, 
which is the generalization of the exponent k 2 /T of the 
previous sections and change the previous assumption 
2ÍIT > k 2 into 2Í7T »wT = r/+/í 2 (l-2r). Therefore we 
get the homogeneous solution A{t) = e^i? -1 ^), (where 
R{t) is the rotation matrix from section lïï"ï]l and thus 



~,out 
X A 
out 

Pa 



B2L I X 



n'" 

Pa 



(A2) 



+e 



-u,t K y/2r r T lx mí^^J o 

+e ' — / dte ^ R(t)[ r ^- , , 

VT Jo '\-Vl - 2rf xL {t) 

where R(T) = 1 was used. The first two lines repre- 
sent atòmic- and light contributions, while the third and 
fourth term account for atòmic noise and light noise re- 
spectively. The light mode, which is naturally mapped 
onto the atòmic sample, is no longer a modulation of the 
upper sideband, as can be seen from the factor (1 — 2r) 
attached to xl(cí, 0) in the second line. Since it is advan- 
tageous to encode the signal at sideband modes, the term 
involving the new disturbed light mode is decomposed 
into a generalization of the familiar plus mode connected 
to the upper sideband 



A as+ 
pin 



,wT 



- 1 



dl 



~R{t) 



p L (ct,0) 
-x L (ct, 0) 

(A3) 



and a small contribution from an orthogonal plus mode 
lying at the lower sideband 



P, 



X 



ls+ 



ls+ _ 



1 



dt e^R{t) 



Ph{ct, 0) 
x L (ct, 0) 



Generalized light modes are denoted by capital letters. 
With this decomposition the atòmic input-output rela- 
tions with noise read 



A 

out 

Pa 



e 2 



az. / x 



A 
r ,in 

Pa 



(A4) 



VT 



,-wT 



k(1 



vin 

pin 
us-\- 



-VT 



,-wT 



1 



TDin 
T 

dte^R{t) 



fxA(t)\ 
fpA(t)J 



-™t K V2r f T , «i . J 



Read-out 

In order to perform the read-out, a second pulse of light 
is sent through the double pass scheme. Subsequently the 
light mode, which is the inverse accented counter-part of 
the mode appearing in the atòmic time evolution (|A2|) 
should be measured. Instead we choose the generalized 
minus mode analogous to the write-in quadratures 1A3|1 . 
The corresponding output quadratures are given by 



Xout 
us- 

T)OUt 
US — 



1 - e 



-wT 



dt 



PL{ct,T) 
-x L (ct,T) 



This can be evaluated by inserting the generalized ex- 
pressions for Pz,(cí, T) and Íl(cí, T). For Pi(cí, T) we 
have 



p L (ct,T) = vT^2f p L (ci,0) + V2rf P L(t) 



x A (t). 



Pl is damped after the first (p-conserving) interaction 
and picks up some noise in return. Subsequently it gets 
some x A -contribution during the second (xlx a -) inter- 
action. The reduced coupling strength k — y/ï — 2r n 
accounts for the damped classical field in the second pas- 
sage. xl(cí, T) on the other hand is given by 

x L {ct,T) = a/Ï - 2f(i L (cí,0) + +V / 27/L(*)· 

xl gets some p A contribution during the first scatter- 
ing interaction i.e. before the relevant trànsits through 
celi walls occur. Subsequently this is damped and ap- 
propriate noise is added. All together both quadratures 
are damped, since both carry the argument (cí,T). This 
means each piece of the pulse contributes after it ran 
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trough the sample twice and has therefore already expe- 
rienced the two relevant trànsits trough celi walls. The 
rest of the calculation is straight forward. In the end re- 
flection losses due to the fourth crossing of a celi wall havc 
to be considered by damping the calculated result by a 
factor y/1 — r and adding appropriate noise terms. The 
resulting input-output relations for the read-out mode 
are 



vout 
us— 

TJOUt 



Cl 



+C 3 



-ce 



Pa 
P 



01 



ls+ 
vin 
A ís+ 

F X A 
F pA 



+c g dtR(t)[e~ 



^US + 

pin 

us-\- 

vir 
pii 

Ui 

FxL 
FpL 

n mt 

e 2 - 



■ c 8 



r>in 
r ls- 

xfc 

FxL 
FpL 



(A5) 



The coefficients c\ to cg can easily be calculated. Since we 
want to focus on the structure of the equation, we don't 
insert these prefactors in order to avoid complicated ex- 
pressions. The new read-out equations differ from (|TT)l by 
the appearance of noise terms (third and fourth line) and 
extra light modes (second line). These contributions are 
small and can be treated as perturbations. ( F x a, F p a ) 
is an atòmic noise mode, while ( F x l , F p l ) , refers to to 
the light mode which is due to the very last reflection. It 
is independent from the light mode ( F p l, F x l ) which 
accounts for the reflections happening between the scat- 
tcring interactions. These intermediate reflections give 
also rise to terms in which only noise associated to xl 
contributes. They are summarized in the expression pre- 
ceded by cg. The appearance of light modes other than 
( X ua+ P us + ) is due to a asymmetry between the plPa 
-interaction and the xlxa present in a realistic setup in 
contrast to the ideal case. The light field has to cross two 
glass walls between the first and the second pass (thus af- 
fecting only the xlxa -interaction). Thus contributions 
from the lower sideband appear and contributions from 
the minus mode do not cancel as in the ideal case. 



APPENDIX B: FIDELITY FOR THE COMPLETE 
STATE TRANSFER INCLUDING NOISE 

In the following subsections the fidelity for storage and 
subsequent retrieval of an unknown state of light will be 
derived for coherent input states and light qubits respec- 
tively. Hereby reflection losses and transverse decoher- 
ence of the atòmic spin state are taken into account as 
explained in section V.A and appendix A. The following 
calculations are based on the input-output relations for 
the complete state transfer. For coherent input states 
they are obtained by inserting equation (|A4|) . which de- 
scribes the atòmic state after a noisy write-in procedure 



into 1|A5(1 . which gives us the final retrieved light state 
in the presence of losses. For light qubit input states 
analogous relations hold. Since we assume an atòmic en- 
scmble at room temperature, no diffusion of the collective 
atòmic mode during the storage time has to be consid- 
ered. Thermal motion of the àtoms ensures, that the col- 
lective mode which was addressed by the write-in beam 
is identical with the collective atòmic mode interacting 
with the read-out pulse during the retrieval procedure. 
Decohering mechanisms such as collisions impairing the 
atòmic state during the storage time oceur on a slow 
seale. As was demonstrated in Q storage times up to 10 
ms can achieved. 



Fidelity for coherent input states 

In order to compute the fidelity for coherent in- 
put states, means and variances of the final quadra- 



tures have to be calculated. 



(X 



fin\ 



and 



can be derived from the expression 
describing the complete state transfer by using the as- 
sumption 2QT 3> wT = r\ + n 2 (l — 2r) (which is a di- 
rect generalization from the approximation 2Í2T 3> k 2 
made in the ideal case) and help of the noise opera- 
tor properties (f x ) = (f p ) = (f x f p + f p f x ) = and 
(/(*)/(*')) = s (t -*')!• The obtained expressions have 
to be inserted into equation (|13|l . which gives the state 
overlap between the input-state to be stored and the final 
state recieved. By considering a gaussian distribution of 
width n for coherent amplitudes the average fidelity can 
be directly calculated as in section lïll C II 



Fidelity for light qubit input states including noise 

The initial qubit state \^i n ) = (a + (3a\ n )\vac) is sub- 
jected to the write-in and read-out procedure which is 
represented by the unitary transformation Un. We ob- 
tain 



finí = U N \V in ) = (a + (3U N al n U f N )U N \vac) 
t 



(a + /3a'f in )UN\vac) 



In contrast to the ideal case, where U\vac) = \vac) could 
be used, Un is a general Bogoliubov transformation. We 
remark that for r — the state transfer can still described 
by a passive transformation. The active contribution is 
cntirely due to reflection losses. This can be understood, 
by noting that reflection losses oceurring between the 
first and the second scattering interaction impair only the 
scattering in the second pass. Therefore the active part 
of the second interaction cannot compensate the active 
part in the first pass as in the ideal case. This leads to 
a term in the generalized atòmic input-output relations, 
which contains only one light quadrature and can there- 
fore not be expressed as a mode-contribution. It plays 
an isolated role in the commutation relations, but adds 
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some extra noise to the variances. Losses due to atòmic 
decay on the other hand are included into the dynamics 
of the scheme in a symmetric way. 

The fidelity for the complete state transfer is given by the 
overlap between the target state |*^) = (a — f3a\ n )\vac) 
and the light state |\£/i n ) which is effectively retrieved 

^ t =|(*^^/™)| 2 =l(^c|(a*-/3*a m )(«+/34jí7 J v|t;a C )| 2 

(Bl) 

a \w, is known, since the input-output relations for the 
complete state transfer are known. They can be written 
in terms of creation and annihilation operators such that 
all occurring modes are independent. The transformation 
is of the type 

n m 

where the coefficients ki, kj are complex numbers. a\ = 
a\ n refers to the state to be stored, while a\ to a^ n rep- 
resent all creation operators which appear in the equa- 
tion, namely contributions from the atòmic input, atòmic 
noise, light-input from the read-out beam and light noise. 
Since we also have noise terms, which cannot be ex- 
pressed as a noise mode (compare equation (|A4|) last 
term) and contributions from the lower sideband (com- 
pare equation (|A5|) ) in which the x- and p quadratures 
are interchanged, we also have annihilation operators in 
this equation which are represented by c\ to cy. Since 
theses contributions are small, they are treated as per- 
turbations to the system. 

The transformation given in l|B2|) can be understood as 
an orthogonal transformation P = P a (g) P c , where P a 
acts on the creation operators and P c acts on the anni- 
hilation operators, followed by an active transformation 



S. With normalization constants N a — l^l 2 an< ^ 

= ^E7=i \k\ 2 ' wh ere N 2 a - N 2 C = 1 and N c < 1, 
(|B2(I can be written as 

n , m 7 

•k-*.(E^-í)+*.(E^) < B3 > 

i=i j=i 
= N a P a a\PÍ + N c P cCl Pl = N a a P + N c c P 
= Vi + N%a P + N c c P = Sa F S* (B4) 

and we have Un = S(P a <S>P c )- In order to compute F qu bit 
from equation l|Bl·|l the expression Un\vüc} has to be de- 
termined. Un\vcic) = S(P a <g> P c )\vac) = S\vac), since P 
is a passive transformation. S on the other hand refers 
to a two mode squeezing operation. As mentioned above 
active contributions are treated perturbatively. The cor- 
responding time evolution S — e N o(o.pc-a p c i ) j g expanded 
in a series to first order and we obtain 

U N \vac) = -j=^=={\ - N c a P cp)\vac) 

By inserting this expression in equation (|B1|I and insert- 
ing the right hand side of (|B2|) for a\ in the fidelity can 
be directly calculated. We find 

In order to obtain the average fidelity we set a = cos(|) 
and (3 — sin(|)e^ and integrate over the whole Bloch 

sphere F qubit = i F qubit {6, 4>) sm(0)d0d(f). The 

results are shown in figures and 
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